UNIVERSITY OF ILLINOIS 
URBANA 



AERONOMY REPORT 





^ (N ASA-CR- 1 42 344) STUDIES OF THE N75- 19889^ 

DIFFERENTIAL ABSORPTION ROCKET EXPERIMENT 
(Illinois Univ.) 96 p HC $4.75 CSCL 04A 

Onclas 

63/46 13466 ) 


STUDIES OF THE DIFFERENTIAL 
ABSORPTION ROCKET EXPERIMENT 


by 


J. C. Ginther 
L. G* Smith 


January 15, 1975 


Library of Congress ISSN 0568-0581 



Supported by 

National Aeronautics and Space Administration 
NGR 14-005-181 


Aeronomy Laboratory 
Department of Electrical Engineering 
University of Illinois 
Urbana, Illinois 




CITATION POLICY 


The material contained in this report is preliminary information cir- 
culated rapidly in the interest of prompt interchange of scientific 
information and may be later revised on publication in accepted 
aeronomic journals. It would therefore be appreciated if persons 
wishing to cite work contained herein would first contact the authors 
to ascertain if the relevant material is part of a paper published or 
in process. 



UILU-ENG 75 2501 


AERONOMY REPORT 
N 0. 64 


STUDIES OF THE DIFFERENTIAL ABSORPTION ROCKET EXPERIMENT 


By 

J. C. Ginther 
L. G. Smith 


January 15, 1975 


Supported By Aeronomy Laboratory 

National Aeronautics Department of Electrical Engineering 

and Space Administration University of Illinois 

Grant NGR 14-005-181 Urbana, Illinois 



11 


ABSTRACT 

Investigations of the ionosphere, in the rocket program of the 
Aeronomy Laboratory, include a propagation experiment, the data from 
which may be analyzed In several modes . This report considers in 
detail the differential absorption experiment. The sources of 
error and limitations of sensitivity are discussed. Methods of 
enhancing the performance of the experiment are described. Some 
changes have been made in the system and the improvement demonstrated. 
Suggestions are made for further development of the experiment. 
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1. INTRODUCTION 

The propagation experiment has been used on all rocket flights con- 
ducted by NASA for the University of Illinois' Aeronomy Laboratory. The 
experiment was designed by the Coordinated Science Laboratory for use in 
the IQSY years of 1964 and 1965. The experiment measures two phenomena, 
differential absorption and Faraday rotation. Faraday rotation is the 
phenomenon of the changing phase characteristics between two propaga- 
ting modes, while differential absorption is characterized by a dif- 
ference in attenuation rates. These two phenomena can be used together 
to determine electron concentrations and electron collision frequencies, 
or independently, by use of suitable approximation, for calculation of 
electron densities. The electron densities used to calibrate probe ex- 
periments, carried by the rockets, are capable of measuring fine 
structure in the electron concentration, but not absolute magnitudes. 

This report is concerned with operation of the differential absorp- 
tion experiment without use of the accompanying Faraday rotation data. 

The objectives are to determine criteria for determining valid data and 

improvement of the accuracy of measurements. It is desired to be able 

-3 

to measure electron concentrations as low as 10 cm . To accomplish this 
the errors in the absorption rates measured on daytime flights must be 
reduced to about .01 dB/sec. For nighttime flights the objective is to 
reduce the errors from values over 1 dB/sec often seen, to values of 
about .1 dB/sec. 

Chapter 2 gives a summary of the theory of the propagation ex- 
periment as it applies to differential absorption. Chapter 3 describes 
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improvements made in the operating and data-reduction systems. Chapter 4 
covers errors in the experiment resulting from unwanted signals which are 
not included in the theory of the experiment’s design. Chapter 5 dis- 
cusses the range of valid data and how this should be defined in future 
work. Finally, Chapter 6 contains suggestions for system improvements 
which can be made without major redesign of the experiment. 
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2. THEORY AND IMPLEMENTATION OF THE PROPAGATION EXPERIMENT 

To gain an understanding of how the propagation experiment is designed 
to operate, the theory developed by Appleton [1932] and Hartree [1931] 
provides an excellent working base. This theory has been improved upon by 
Sen and Wyller [1960] but the Appleton-Hartree equations provide the 
necessary tools to understand the experiment and are the basis upon which 
the experiment was designed. 

The following sections will present an outline of the Appleton-Hartree 
theory and how this is changed by Sen and Wyllex. Then a description of 
how the experiment is designed and an outline of the data reduction pro- 
cess will be presented. 

2.1 The Theory of Appleton-Hartree Under Quasi-Longitudinal Propagating 

Conditions 

In a region of free electrons, such as the earth's ionosphere, where 
there is a magnetic field present the medium becomes bi-refringent. This 
means that a traveling wave entering this medium will be split into two 
elliptically polarized modes, each mode subject to a different index of 
refraction. The theory of Appleton-Hartree determines the modes of prop- 
agation and the indexes of refraction of these modes in the earth's 
ionosphere under the following assumptions: electron collisions with neu- 

trals are independent of electron energies, the medium of propagation is 
electrically neutral with a uniform charge distribution, the magnetic 
field is uniform throughout the medium, and the ions, because their masses 
are much greater than that of an electron, are stationary. The theory 


follows: 
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Assume a vertically propagating wave which has the form 

= E q exp [ j ( tot - kz)] (2.1) 

k is the propagation constant, w is the operation frequency, and E a 
reference field 


E 

o 


- E x + E 
x y 


y « 


( 2 . 2 ) 


The polarization of this wave is defined as 



(2.3) 


If the fields are in the form of equation (2.2) they will have 
refraction indices given by 


n 


= 1 + 


x 


E E 
o x 


= 1 + 


E E 
o y 


wJ~E 


0^0 


(2.4) 


In the above, P and P are the polarization fields in the x and y direc- 
^ if 

tions. The motion of the free electrons in the ionosphere will be governed 
by the equation of motion: 


ra* + mr_ = - e (r x 5) - eE_ . (2.5) 

v is the electron-neutral collision frequency, and B the magnetic flux. 

The polarization can be found from the electron displacements as 



where N is the electron density. 


Before continuing the development of the theory, the following table 
of standard symbols is presented. 





is the plasma frequency and is defined as 



is the gyrofrequency defined as 


V-H 



which is divided into two conponents, one along the direction of propa- 
gation and one transverse to the propagating wave. 

With the above set of symbols, equations (2.5) and (2.6) can be com- 


bined to form 



6 


where 


E E = P 
o — — 


( 2 . 7 ) 


M' 1 = - X~ l 


\ 


1 - jz 



1 - 3 m Z 

JX T 

0 


~ tjX n 


1 - jz 


When this result is combined with Maxwell's equations, the characteristic 
equation for the polarization becomes 


Y 2/y 

9 i /7l / i J' 

* 1 = 0 

which when solved for the polarization yields 


R = -3 


y t /2Y l - L J / / 4 I I 2 

+ /I + 


-X-3Z ' yf* ' 1+X-jZ 

and the solution for the refractive index is 


( 2 . 8 ) 


( 2 . 9 ) 


n 


2 


1 - 


X 

1 - jz - jY l r 


( 2 . 10 ) 


which is known as the Applet on -Hartree formula. 

For rocket flights from Wallops Island, with which this study has 
been primarily concerned, the equations can be simplified since the direc- 
tion of propagation is approximately along a magnetic field line so that 
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V 



<< 


{ 

1 - X - jz 


is valid. 

When this approximation is used inequations (2.9) and (2.10), the 
results are 


R = ±3 


( 2 . 11 ) 



1 - 


X 

WZ ± i L 


( 2 . 12 ) 


The solution for R indicates two circular modes will be propagated. The 
+«7 is known as the extraordinary mode and the -q as the ordinary. In the 
refractive index relation, the plus sign on is associated with the 
ordinary mode and the minus with the extraordinary mode. 

A further approximation in the equations may be made on the basis 
that the plasma frequency in the D and E regions where the experiment is 
intended to be used, is much smaller than the frequency of the propagating 
wave used in the experiments, thus 

X « 1 (2.13) 


n - 1 



1-j'Z + x L 


.(2.14) 


Using (2.13) in (2.12) yields 
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The index of refraction is divided into its real and imaginary parts 


n = m + j x 


n = 11 - 


U (1 ± * ) 


- Jl 


\xz 


.Cl ± * £ ) 2 + z 2 


(2.15) 


and the equation for the propagating wave may be written as 


S_ ■ exp(-x^Z) exptjCwt-yk^Z] 


(2.16) 


This is the equation of a wave which decays in amplitude as it propagates. 
It is this phenomena described by the coefficient of absorption 


K s \ X (2.17) 

with which this study is concerned. 

2.2 The Sen-Wyller Theory 

An improvement on the Appleton-Hartree theory was developed by Sen 
and Wyller [I960] to correct the assumption that the collision frequency 
is independent of electron energy. The collision frequency is assumed to 
have the form 


v = /(d) (2.18) 

/(D) is a function of the electron velocity which is assumed to be 
Maxwellian. The theory is quite complicated and will not be reproduced 


here. 



9 


The result of the. theory is to give an index of refraction 


n 


2 2 4 2 2 

4 + Bsinii±B sin i|j - C cos ih 

. 2 

D + E sin ^ 


C2.I9) 


^ is the angle he twe enr _ th e' wave^n orma 1 and B, with the other variables 
defined as follows: 


A = 2 C J (C I + t ni ) 


B + t'lll ^ 


C - 2 


D = 2 l 


E = 2 l 


III 


tj, and £jjj are the dielectric tensor elements defined as 


tj. = (1 -a) -3 b 


l II = + Ho-e)/2 


l III = C^-C^+eD/2] + 3 [b - C/+4)/2] , 


in 


which a, b, c, d, e, and f are functions of the 0 script integrals 
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Figure 2.1 shows how the Sen-Wyller and Appleton-Hartree theories 
differ, especially at high collision frequencies found at the lower alti- 
tudes. It is clear from this that the full Sen-Wyller theory will be 
necessary for analysis below 85 km. 

2.3 Operation of the Experiment 

In order to measure the absorption properties of the atmosphere, the 
two modes are transmitted from the ground to a rocket as it passes through 




h (km) 


Figure 2.1 Appleton-Hartree versus Sen-Wyller differential absorption index per 

unit electron concentration as a function of electron collision frequency. 
The approximate Appleton-Hartree curve is normalized to agree with Sen- 
Wyller curve for vanishing collision frequency [Meohtty et 1967]. 


/ 
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the ionosphere. By measuring the differences in the rate of absorption of 
the two modes (differential absorption) and combining the results with the 
Faraday rotation experiment, or by using a collision frequency model, the 
electron concentration of the ionosphere can be measured. 

A block diagram of the original system is shown in Figure 2.2. The 
signal generated by the configuration consists of an ordinary wave at fre- 
quency f - 6, where 6 is 250 Hz, and an extraordinary wave at f o + 6. 

The resultant wave can be visualized as an electric vector of frequency / 
spinning in the ordinary direction with a frequency / - 6. The magni- 
tude of the vector is governed by an ellipse spinning at frequency <5 in 
the ordinary direction as shown in Figure 2.3. The above assumes the 
ordinary magnitude is larger than the extraordinary and that <5 is 
positive. 

Analytically, the result can be shown by taking the ordinary wave to 
be 

A A 

i Q sin[(w-6 ) t - k^Z] - j Q cos[(w-S)i - 
The extraordinary wave is then taken as 


i R sin[(u)+<5)t - k Z] * q R cos[(w+S)t - k Z] 

X X 

k and k are the propagation constants for the o (ordinary) wave and the 
x (extraordinary) wave. 

The unit vector for the rocket antenna spinning at may be written 
as 

A A 

a = - i sinw t + 3 costu t 

— V V 


( 2 . 21 ) 
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Figure 2.2 Original system block diagram [Knoebel et al. 3 1965]. 
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Figure 2.3 Generation of polarization ellipse (adapted from Salah and Bowhill [1966]). 



Now by substituting in Z = being the rocket velocity, the total 

signal received by the antenna is proportional to 


E • a = ' ( 2.22 ) 

-[Q sin(w-<S-& v )t + R sin(cu+<5-fc v )t] sinw t + \-Q cos (w-6-fc v )t 
v Of xv v Of 


+ 1? cos (w+S-fe V )t] costo t + -Q cos(u-6-fc V -u> )t + R cos(w+S-k V +x )t 
Xf f Off x f f 


Making the substitutions 


A = (-6 - k V - a) )t 

Off 


= (6 - k v + 
x f 


uj )i 

ji ' 


The output can be written as 


-Q[cosA cosajt - sinA sinujt] + i?[cosB coswt - sinF sinwt] 


= [ -Q cosA + R cosB] cos xt + [4? sinA - R sinB] sinwt 


(2.23) 


By employing e 


Obit 


= cosoj t + o sinwt the above can be written 


Re[{-Q cosA + R cos B) e^'^] + Im[(Q sinA - R sinF) e t?au ’] 
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The receiver detects the amplitude of the coefficient of e which is 

2 2 1/2 
\QT * R - 2 QR(A - B)] 

which may be approximated by 

1/2 

(Q 2 * R 2 ) [1 - -Af cos (26 - (k - k )v + u )t] 

Q^+R* 

The resulting output can be seen to consist of a dc component and one at 
26 - (k - k )v + 2to rad/sec. 

The receiver output is telemetered to the ground and transferred to 
the transmitting van over data lines. At the van, the dc and 26 + 2u^ 
components ((-k + k^)u^ is assumed small) are compared in the percent 

modulation detector. The modulation is kept at a constant percent, and 
thus the axial ratio of the ellipse at the rocket is kept constant, by 
adjustment of the transmitted powers of the extraordinary and ordinary 
waves. This is accomplished by a servo motor which is controlled by the 
output of the modulation detector. The servo moves a coil within a wave- 
guide beyond cutoff. This acts as an attenuator which is logarithmic in 
voltage and thus linear in dB. The position of the potentiometer, which 
is used to measure the attenuation factor, is then recorded to give a 
record of the rate at which power was added to the X wave. The two 
methods of getting the potentiometer information to the recorders, AM 
and FM, are discussed in Chapter 3. Figures 2.4 and 2.5 show examples of 
the outputs on a chart record. 




LBE MK 65 


55 

TIME FROM LAUNCH (sec) 


Figure 2.4 Section of chart record for Nike Apache 14.513 showing differential absorption 
signals in relation to lower band edge (LBE) , band center (BC) and upper band 
edge (UBE) . This illustrates the AM system. 



Figure 2.5 Section of chart record for Nike Apache 14.511 showing differential absorption 
signals in relation to lower band edge (LBE) , band center (BC) and upper band 
edge (UBE) . This illustrates the FM system. 
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2.4 Procedures for Recovering Rates of Differential Absorption 

Recovery of rates of absorption data to the time of this writing, 
has made use of CDC 1700 and IBM 360 computers. The 1700 is being re- 
placed and its portion of the reduction process will also be circumvented 
by use of digital tapes prepared by NASA, since the 1700 is not used in 
actual reduction but as a means of getting data. from the analog to the 
digital form. 

The first step in the process of rate recovery is to convert the 
information contained in the analog tapes to a digital format suitable 
for. analysis on the IBM 360 computer. Figure 2.6 is a block diagram of 
the system used for this process. The signals representing the time of 
data, position of the extraordinary power potentiometer and a 100 kHz 
reference signal are desired from the tape. The 100 kHz signal is used 
to measure any tape speed variations and feeds a corrective signal to 
the discriminator processing the extraordinary power signal. Errors in 
amplitude in the time signal do not cause errors in the time interpreted, 
so there is no need to correct this signal. After detection, the time 
signal is still in a modulated form and is processed by the time code 
demodulator which sends the proper signal to the digitizing unit. 

The output of the extraordinary discriminator is either sent directly 
to the digitizer in the cases where the FM system of potentiometer posi- 
tion recording is used or to a rectifier when the AM system is employed. 
The rectifier detects the amplitude and this is then transferred to the 
digitizer. The outputs of the digitizer are placed onto digital tape by 
the CDC 1700 computer. This tape is used with the IBM 360 for all further 
processing. 




Figure 2.6 Block diagram of digitizing process. 
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On the 360, DACAL is run first. It processes the data from the pre- 
flight calibrations to determine which signal levels correspond to given 
attenuator settings. 

The final step of recovering the rates of absorption is accom- 
plished through program DAPROC. This program analyzes the data on the 
tape and gives an output of rates of differential absorption in dB/sec. 
The program may also be directed to punch cards with the differential 
absorption rates and complete trajectory data for the rocket. This is 
needed as input for the program which analyzes the data to determine 
electron concentrations. 

2.4.1 Computer determination of electron concentration from dif- 
ferential absorption rates. Mechtly et at. [1967] discuss the evalua- 
tion of electron concentrations when Faraday rotation data are available. 

We will deal here with the evaluation of these data in the absence of 
Faraday rotation. 

If the rocket trajectory and the characteristic of the magnetic 
field are known the only unknown in the equations for determining how 
each mode is absorbed are collision frequency and electron concentration. 
To solve for the electron concentration, a model for the collision fre- 
quency must be used. It is assumed that the collision frequency is pro- 
portional to pressure. 

\ = C x 10 5 pO) (2.25) 

The problem of determining C is discussed by Mechtly [1974] . 
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With the value of C chosen, and given pressures taken from 
1972, the necessary information to solve for the electron density is 
at hand. The program assumes an initial electron density and then 
applies the exact form of the Sen-Wyller equations to determine a dif- 
ferential absorption rate. This is compared with the measured rate and 
the assumed electron density modified. This process is repeated six 
times after which it usually has converged to within better than one 
percent. 
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3. IMPROVEMENTS MADE IN DATA ACQUISITION AND REDUCTION SYSTEMS 

As stated in the Introduction, one of the objectives of this project 
was to find ways to improve the data, specifically at low altitudes. In 
an effort to improve the operation of the data acquisition and reduction 
systems, several of the components making up the differential absorption 
experiment have been improved. This chapter will examine why the AM data 
link was replaced by the new FM system and the resulting improvement. 

Then the operation of tape speed compensation will be reviewed and its 
effects on the data discussed. Finally the possible use of filtering to 
recover data from the AM system will be reviewed. 

3.1 FM Data System 

3.1.1 Theoretical superiority of frequency over amplitude modula- 
tion as a data link. Part of the original equipment in the University of 
Illinois van consisted of an amplitude modulated link to transmit the 
information about the attenuator position to the recording station. 

This system, shown in Figure 3.1, simply generated a signal with an 
amplitude proportional to the attenuator position. Due to age, the per- 
formance of the four oscillators was degraded and replacements became 
necessary. With the necessity of replacement, the possible use of a fre- 
quency modulated system was examined to see if improvement could be 
made . 

For transmission systems, the figure of merit may be used as a; 
basis of comparison. The figure of merit is defined as the signal-to- 
noise ratio out of a receiver, divided by the signal-to-noise ratio of 
the incoming signal. 




Figure 5.1 Block 
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Y 


S fN 
o o 

s./n u 

v M 


(3.1) 


with N m being the noise power within the receiver bandwidth. 

For an AM system with a signal Am(t) cos 2 vf t, y was found by Taub 
and Schilling [1971] to be 


Y = 


m 2 (t) 


1 + m it) 


(3.2) 


2 

where m (t) is the mean square value of the modulation coefficient. 

For the period of the flight when the differential absorption is 
just beginning, the signal can be approximated by a sinusoid of ampli- 


tude C leading to a figure of merit 


m ' 2 * C 2 


(3.3) 


For an FM system under the influence of a sinusoid, Taub and 
Schilling show that 


‘ T PW -!* 2 < 5 - 4 > 

B is the modulation index for the FM signal, f/f M > f m being the maximum 
frequency deviation. Since the same potentiometers would be used in 
the FM system it will have the same modulation index as the AM system. 


B = C. With this 
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(3.5) 


To compare the two systems the ratio of their figures of merit will be 
used 


3 

]FM = _2_ 
y AM 


This indicated that the implementation of an FM system will improve the 
noise rejection of the data transfer process by more than three. 

Because of this result, an FM link was installed. 

3. .1.2 Design of frequency modulated data link. Figure 3.2 shows 
how the FM system was implemented. To minimize alterations to the sys- 
tem, as few changes as possible were introduced. In place of the four 
audio oscillators of the AM system, a single precision dc supply was 
installed. The voltages present at the variable contacts of the poten- 
tiometers are then used as the inputs to voltage controlled oscillators. 

To choose the VCO frequencies, tests were run on the data lines at 
Wallops Island, Virginia, to determine their frequency response charac- 
teristics . Figure 3. 3 gives the frequency responses of four data lines 
and the locations of IRIG proportional bandwidth channels 1 through 10. 
It is desirable to use the highest useful frequencies since the data 
rate increases with the higher channels. Under these restrictions, 
channels 4 through 7 were selected. 


6 ^ 


= f(2 


C 2 ) = 


3 + 


3 fZ 

y C 


(3.6) 




Figure 3.2 Block diagram of frequency modulated data link. 
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3.1.3 Data reduction with frequency modulated system. As noted in 
Chapter 2, with the addition of the FM system, the process of digitiza- 
tion is simplified since the additional rectifier in Figure 2.6 is no 
longer required. The change from the AM to the FM system, however, 
revealed an inefficiency in the median value subroutine of DAPROC. 

When 14.511, the first rocket launch during which the FM system was 
used, was processed it took DAPROC four minutes of IBM 360 computer time 
to analyze 4 seconds of data, where previously up to 40 seconds had been 
processed when the AM system was employed. To show how this happened. 
Figures 3. 4, .and 3.5 shojw how samples of data for both the AM and FM sys- 
tem, might look. The noise present in the AM system is a function of 
two factors:, 1) increase noise in data transfer, and 2) incomplete 
filtering of the carrier signal by the rectifier. 

If the signal has no noise and the curve is monotonically increas- 
ing, the median value will be located at the mid-point on the time axis 
of the period of interest. With the noisy data, the routine would find 
the median in a reasonable time because a point of magnitude equal to 
the median is found early in the data. If the data are smooth, as that 
generated by the FM system, the routine will search through nearly 
half the data to get to the median. With 10,000 data points, this means 
5,000 passes through the data making comparisons, or 5 x 10^ comparisons. 

To rectify the problem, the statement directly above 200 in DAPROC, 
Figure 3.6, is changed to Test = XPWR(ISTOP/2) , moving the initial 
search point to the mid-point of the time slot. With this change, 

DAPROC was able to process 40 seconds of data from 14.511 in less than 


three minutes. 



SERVO POSITION 


30 



TIME 


Figure 3.4 Possible AM system data. 





SERVO POSITION 



TIME 


Figure 3.5 Possible FM system data. 




START 


II STOP = 
I STOP+1 


NN2 = 

I STOP/ 2 


MN= 1 


TEST « 
XPWP(l) 


NET *0 
NGT =0 
NEQ °0 
K * MN 


NGT? 
NGT +1 


NEQ 8 
NEQ +1 


XPGT4NGT) 
*XPWR( K ) 


XPLT ( NLT) 
= XPWR (K) 


K = K+i 


JJ STOP, 


Figure 3.6 Flowchart of original DAPROC program. This flow 
chart is a corrected version of that given by 
Stekys and Mechtly [1970] . The program listing 
in that reference is correct and is currently 
stored on the IBM 360. 
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Figure 3.6 Continued 
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3.1.4 Verification of data link improvement by use of FM system. As 
a means of testing the improvement made by using the FM system, data, 
digitized while the potentiometers at the University of Illinois van were 
locked, were analyzed as if it represented real data. Table 3.1 lists the 
rates of differential absorption that resulted from noise in the systems, 
(During this test tape speed compensation discussed later was employed.) 
The standard error of the FM system is ,006 dB/sec with a maximum change 
of .01 dB, the minimum change the programs are designed to measure. The 
older AM system showed a larger .037 dB/sec standard with errors up to 
.09 dB/sec possible. 

As a result of the FM system, the data transfer can be eliminated 
as a factor in the overall noise of the measurements involved in the 
differential absorption experiment. 

3.2 Effects of Tape Speed Variation on Data and its Correction 

Figure 3.7 shows a chart record of the output of a discriminator 
after being passed through a low pass filter of 10 Hz to remove the 
real data, which are on a carrier at 145 Hz. What this signal repre- 
sents is tape speed noise caused by flutter in the tape playback system. 

In the analysis of data one second of data at a time has been used in 
most cases. If the tape speed noise is assumed periodic, the worst case 
occurs when one second of data corresponds to a multiple of the period 
of the noise plus a half period. 


1 sec = NT + T/2 
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Test 

TABLE 3.1 

of Data Transfer System 

Referenced 

D.A. Rate in 

dB/sec 

Time 

AM system 

FM system 

1 

-.02 

0 

2 

-.09 

1 

• 

o' 

H 4 

3 

-.04 

0 

4 

0 

.01 

5 

.02 

0 

6 

0 

0 

7 

/ .02 

0 

8 

0 

0 

9 

0 

0 

10 

-.04 

0 

11 

.07 

0 

12 

.02 

0 

13 

-.02 

-.01 

14 

0 

.01 

15 


0 

16 


.01 

17 


-.01 
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Under this condition, the total error caused by a signal of frequency / 
and amplitude A is 


Error 



in the one second of data. Since two adjacent seconds are used to com- 
pute differential absorption rates and their errors will be in opposite 
directions, the error in measurement will amount to 


Error = —4 . 

The noise in Figure 3.7 was found to have a primary frequency of 
about 3.5 Hz and amplitude of .4 dB, so the resulting error would be 
.02 dB. 

3.2.1 Implementation of tape speed compensation . In order to 
compensate for variation in the speed of tape passing through the play- 
back system, tape speed compensation was implemented. TSC operates 
by recording a 100 KHz signal onto the data tape at the same time the 
data are recorded. During playback, this signal is fed into a reference 

discriminator centered at 100 KHz. The output of the discriminator is 
proportional to the difference between 100 KHz and the frequency of the 
signal from the tape. 

The signal to be digitized is passed through a delay line and then 
into the discriminator. The delay line compensates for the time it takes 
the 100 KHz signal to be processed in the reference discriminator. The 
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normal output of the discriminator is then added to the signal from the 
reference. The resulting signal will have noise due to the tape speed 
variation reduced by 30 dB. 

3.2.2 Effectiveness of tape speed compensation in reducing noise. 

To see what kind of improvement tape speed compensation would make, 
flight 14.514 was reprocessed both with and without TSC. Table 3.2 
shows the results of these runs. The standard correction is ±.015 dB/ 
sec. For a correction of this order, the 30 dB noise reduction of the 
TSC system is sufficient to reduce tape speed noise below the levels 
measurable by the programs. 

3.3 Filtering of Signals to be Processed from AM System to Improve Data 

Another test performed, to improve the quality of data processed 
for flights employing the AM system, involved the use of filtering. A 
bandpass filter is placed between the discriminator output and the 
rectifier in Figure 2.6. The filter will remove noise that is not 
associated with the data signal, such as dc offset noise. The results 
of this test for flight 14.361 are shown in Table 3.3. The correction 
has a standard value of .02 dB/sec. This is a small error relative to 
those seen in the data, and thus does not warrant reprocessing of 
older flights. 

3.4 Results of Improvements 

All of the new techniques discussed in this chapter have contributed 
in some way to the more exact evaluation of data, but none accounts for 
the large discrepancies noted in Chapter 1. Possible explanations for 
the larger errors are covered in the following chapter, along with 
errors that do not affect t data at low altitudes but become important 
higher up in the ionosphere. 



TABLE 3.2 

Test of Tape Speed Compensation on Data 
from 14.514 
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Differential Absorption Rates (dB/sec) 


TIME 


no 

Difference 

after launch 

TSC 

TSC 

TSC - no TSC 

40 

0 

0 

0 

41 

0 

.02 

.02 

42 

-.11 

-.11 

0 

43 

.08 

.07 

-.01 

44 

.05 

.05 

0 

45 

0 

.02 

'.02 

46 

-.15 

-.16 

-.01 

47 

.10 

.10 

0 

48 

.13 

.13 

0 

49 

-.05 

-.03 

.02 

50 

.11 

.11 

0 

51 

.02 

0 

-.02 

52 

.02 

.03 

.01 

53 

0 

0 

0 

54 

.59 

" .57 

-.02 

55 

.16 

.16 

0 

56 

.72 

.70 

-.02 

57 

1.12 

1.16 

-.04 

58 

1,25 

1.25 

0 

59 

2.34 

2.33 

-.01 

60 

3.09 

3.09 

0 

61 

2.86 

2.88 

.02 

62 

2.67 

2.68 

.01 

63 

2.49 

2.53 

.04 

64 

2.23 

2.22 

-.01 

65 

1.36 

1.36 

0 

66 

1.42 

1.41 

-.01 

67 

1.64 

1.65 

.01 

68 

1.45 

1.44 

-.01 

69 

1.13 

1.12 

-.01 

70 

-.02 

-.03 

-.01 

71 

-.08 

-.06 

.02 

72 

.32 

.29 

-.03 

73 

-.11 

-.01 

.02 

74 

.94 

.92 

-.02 

75 

1.02 

1.01 

-.01 
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TABLE 3.3 

Filter Test on 14.361 


Differential Absorption Rates 


Time 

after launch 

No Filter 

(dB/sec) 

Filter 

Change 

48 

.16 

.20 

-.04 

49 

.48 

.46 

.02 

50 

-.11 

-.11 

0 

51 

-.36 

-.37 

.01 

52 

.40 

.42 

-.02 

53 

.39 

.37 

.02 

54 

-.28 

-.28 

0 



0 
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4. ERRORS IN THE SIGNAL RECEIVED BY THE ROCKET 

This chapter will deal with errors not connected with the processing 
or recording of data, but with errors at the rocket itself. The first 
two sections present a more detailed analysis of the electric field than 
has previously been undertaken and shows some faults in the experimental 
assumptions. The last section deals with actual errors in the receiving 
system itself. 

4.1 Theory of Polarization Errors 

The type of signals which can be expected from errors in the pol- 
arization, those due to imperfect generation of the extraordinary and 
ordinary modes, are examined. The wave form detected by the receiver 
was derived as follows. 

Two signal components are considered to be present at the frequency 
to - 6 . They are the ordinary wave , 

i Q sin[(w-S)£ - X Z] - j Q cos[(w-fi)« - K QQ Z] 
and the unwanted extraordinary wave, 

i S sin[(io-6)i - X^Z] + j cos [0-6)* - 

At the frequency to + 6 there are also two components present. The 
extraordinary wave is represented as 

i R sin [0+6)* - K xx Z] + J R cos[0+<$H - K^Z] 
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and an unwanted ordinary wave is 


£ P sin[(a>+5)£ - K ] - j P cos[(<o+6)£ - X^Z] 


The unit vector for the antenna is taken as in Chapter 2 to be 


a = -£ sinco^f + j cosu^t (4.1) 


and the displacement is taken as 


2 = (4.2) 

X , X , X , and X represent the propagation constants. The first 
subscript represents the mode of the wave, o for ordinary and x for 
extraordinary. The second subscript is the frequency of operation, o 
for the correct ordinary- frequency (u>-6), and a: for the correct extra- 
ordinary frequency (to+6) . v p is again the rocket velocity and the 
spin rate of the rocket. 

a 

The signal received, represented as E • a is 


-[Q sin (co-6-X^ypt + S sin( u-S-K^v^t 


+ R sinCdj+S-X^upt + P sinCw+S-X^v^) t] • s intuit 
+ [-Q cos (a + S cos( w-6-K^v^t 
+ R cos(w+6-X y^t - P cos (<o+6 -X^y^,) t] cosco^t 
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which is 


- Q cos(co-cS ~K y -to )t + S cos y -to )t 
oo r v xo r r 


+ R cos(to+6-^ y +to )t - P cos (to +6 -if y -to )t 
xx r v ox v r 


This may be rewritten as 


[- Q cos (-6-# y -w )t + S cos (-&-K y +<o )t + P cos (S-K v +<o )t 
oo r r xo r r J v xx r v 


- P cos (6-P y -to It] cos tot + [Q sin(-6-.K y -to )t - 5 sin(-6-P V +<o )t 
ox r v J J L oo r v v xo v v 

- R sin(6-P y +(o )t + P sin(<5-Z V -to )tl • sintot 
v xx r r J - ox r v } 1 

or in complex notation as shown below 

Re{[-Q cos (-S-X v -io )t + S cos(-6-P y + to )t 
oo v v v xo r r J 

+ R cos (6-.K y +io )t - P cos(6-K y -to )t 
xx r r K ox r r 

-jQ sinC-S-^y^-to^t + jS sin 
+ SR sinCS-^y^+to^t - j‘P sinCS-Z^-u^) t]e' 7 "“*} 


The amplitude detected by the receiver is 
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{Q 2 + S 2 + R 2 + P 2 - 2QS cos[-(Z -K )v - 2ai It 

1 ^ oo xo J r r J 

-2QR cos [-2 6-(JC -*)» - 2io Jt + 2QP cos [-26 -(K -K )y It 

oo xx v r* J 1 oo ox J 

(4.3) 

+2SR cos [-26 - (■^ xo ~^ xx ) v p ] t + 25P cos [-26- + 2^]t 

1/2 

-2i?P cos[(Z -Z )y + 2oi It} 

LV xx ox J r r J 

The relation 


(1 + X) =1 + | 

for X small may be used to reduce expression (4.3). Table 4.1 lists the 
resulting frequencies and their proportional amplitudes. 

Each frequency can be given a physical interpretation to help under- 
stand what errors are occurring. Frequency b at 26 + 2w is the fre- 
quency discussed in Chapter 2. Frequencies a and f are caused by the 
error signals at each transmitted frequency combining with the correct 
signal to form a small linear signal. The linear signal’s received 
strength is then dependent on the antenna orientation and is therefore 
modulated at 2^. Frequencies o and d are both caused by signals of 
the same polarization, but different transmitter frequencies. Since 
they are the addition of two circular modes in the same direction, they 
result in a circular mode and are unaffected by rocket spin. The 26 is 
caused by the envelope generated by the difference in frequency. 
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TABLE 4.1 



Received Frequencies 

Including Polarization 

Errors 

a 

(Z -Z ) + 2lo 
v oo xo J r 

2<o 

r 

QS 

b 

26 + (Z -K )v + 2u 
v <%» rear r j 5 

26 + 2co 

r 

QR 

c 

26 + (X -Z )y 
K oo ox J r 

26 

QP 

d 

26 + (Z -Z )y 
a:c> xx T 

26 

SE 

e 

26 + (Z -Z )i> - 2<o 

v xo ox J v r 

26 - 2 w 

T 

SP 

£ 

(Z -Z )y + 2w 
v xx ox J r r 


KP 


dc 


[P 2 +Q 2 +R 2 +S 2 ] 


1/2 
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Frequency e is generated as in Chapter 2 but since the ordinary wave is 
now higher in frequency, the ellipse rotates opposite to the direction 
in Chapter 2 and the rocket spin thus subtracts from the total. The 
dc signal is the average signal strength present at the rocket. 

4.2 Effect of Polarization Errors on System Response 

Figures 4.1, 4.2, and 4.3 show the probe current and DA determined 
electron concentrations. Note that in all cases the DA data from the 
2.225 MHz system indicate low electron densities above 80 km. To see 
if this could be due to polarization errors, a Fourier transform pro- 
gram developed by K. L. Miller [Eduards , 1973] was used to recover the 
power in the frequency components at 26, 26 + 2 and 2io y . These powers 
are shown in Figures 4.4 to 4.6. To see what has happened, it is 
assumed that the error extraordinary signal is negligible; thus, S = 0 
in Table 4.1. The amplitude of the ordinary error signal, P, may not be 
assumed to be 0 since as the extraordinary wave is absorbed the power to 
the extraordinary transmitter is increased, and thus the power in the 
error wave, P. As the ratio of ordinary to extraordinary is being kept 
constant, the magnitude of the error ordinary signal continues to in- 
crease. At some point the magnitude of the error ordinary wave and the 
extraordinary wave at the rocket are equal. At this point, the servo 
system will follow the 26 rad/sec signal generated by the two ordinary 
modes rather than the correct 26 + 2w r rad/sec signal. Appreciable de- 
gradation of data is noted to occur when the 26 rad/sec signal power 
reaches about one-fifth the magnitude of the 26 + 2ui^ rad/sec component. 
To get a measure of how much useful data can be extracted from a flight, 
the powers of the 26 and 26 + 2id^ components at low altitudes can be 



ALTITUDE (km) 


PROBE CURRENT (A) 


NIKE APACHE 14.511 

PROBE CURRENT 

a 2.225 MHz DA. 
a 3.385 MHz D A. 

C - 7.3 x I0 5 



ELECTRON DENSITY (cm -3 ) 


Profiles of probe current (upper scale) and electron density (lower 
scale). Data from the differential absorption experiment alone 
are shown. 
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io io 2 


ELEC 

Figure 4.3 Profiles of prc 
Data from the c 



WER (arb. units) 










compared. The difference in power in decibles minus seven decibles gives 
the total range of useful absorption, the sum of the rates for each 
second. Table 4.2 lists the power differences and useful data range 
for 14.511, 14.513, and 14.514 and the second of flight and altitude 
where the limit is passed. 

4,3 Reflected Waves 

The following will give an idea of the error to be expected when 
reflections are present. The reflection coefficients for the extra- 
ordinary and ordinary waves, as measured from the rocket, may be repre- 
sented as T and T . Since the reflection coefficients are measured 
x o 

from the rocket, they include absorption over the path of the reflected 
wave. The propagation will be along a magnetic field line so that the 
mode of the reflected wave is the same as that of the incident wave. 

The error in polarization is assumed to be zero. 

The total electric field may be calculated as the sum of four 
propagation vectors. First there are the two upward propagating modes 


i,Q sin[(u>-6)£ - KZ~[ - jQ cos[(u)-<5)£ - KZ\ 

and 

iR sin[(w+5)t - KZ] + jR cos[(io-6)t - KZ ] 

where it is now assumed for simplicity that the propagation constant, 

K, is the same for both modes. The two reflected vectors can be repre- 


sented as 
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TABLE 4.2 

Power Differences, Useful Data Range, and Limit 
of Data in Time and Altitude for Three Rocket Flights 


Flight 

Difference 
Power (dB) 

Useful 
Range (dB) 

Last 
Time (sec) 

Useful Data 
Altitude 

14.511 

20 

13 

59 

75 

14.514 

19 

12 

61 

75 

14.513 

19 

11 

59 

75 



r>s 


and iRY sin[(w+6)t + KZ] + jRT cos [(or 6)* + KZ]. 

O J/ 

By going through the procedure of Section 4.1 the amplitude of the 
signal detected by the receiver is 


{[1 + r o 2 ]Q 2 + [1 + t x 2 1 r2 - cos (26 + 2 

+ 2Q 2 T cos (2 Kv )t - 2QRT cos [26 + 2 Kv + 2w )i 

0 kC ^ 


- 2QRV cos (26 - 2Kv + 2 ujt 

O V V 

2 1/2 
♦ 2irr cos(2K« - 2«j?r r cosC26 + 2» )t> 

# v r (9 0? T 

If r and r are small and R is much less than Q the expansion used 
x o 

in Section 4.1 is valid and the amplitude and frequencies present in the 
output are shown in Table 4.3. Thus, under the conditions of reflected 
waves with absorption the modulation index, M , is 


QR( 1 + 


M = 


r r ) 

o x J 


[d + r 2 )Q 2 * (l + t x 2 )r 2 ] 1/2 


(4.4) 


It can be seen that changes in the shape of the ionosphere above the 
rocket will cause the index M to change and the system to wander looking 
for a proper setting to hold M constant, as shown in Figure 4.7. 

When perfect reflections occur the expansion used breaks down. To 
examine this case let R = 0, so that only the ordinary wave is present, 
and let r = -1. The receiver output is then 
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TABLE 4.3 

Frequencies in Receiver Output when Reflections are Present 

for T r « 1 and R « Q 
x> o 


Frequency 

dc 

2Kv r 

26 + 2to 

r 

26 + 2(o + 2 Kv 

r r 

26 + 2 to - 2 Kv 
r v 


Amplitude 

22 22 X/2 

[Ci+r + (1+r.pirj 

2 2 

r $ + r F 

O X 

qr + rrffl 

O X 

r *« s 

r o e» 




TIME FROM LAUNCH (sec) 

Figure 4.7 Section of chart record for Nike Apache 14.520 showing differential absorption 
signals in relation to lower band edge (LBE) , band center (BC) and upper band 
edge (UBE) . This is the FM system used for a nighttime flight. 


On 
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[2 Q 2 - 2 Q 2 cos (2 Kv r t)] 1/2 


This can be rewritten as 


2Q\ sin(Ku^t) | 

The frequencies in the output are shown in Table 4.4. Note that the dc 
terra is larger than that for only an upward traveling ordinary wave when 
it would be simply Q. 

It should be noted that this analysis does not hold near the ordinary 
reflection point the wave is no longer circular, but linear. Also the 
reflection occurs at X = 1 due to ray bending. 

To connect the results of sections one to three of this chapter 
Figure 4.8 shows the output spectrum of the receiver during flight 14.270. 
All of the frequency components around 500 Hz predicted are present. 

4.4 Receivers 

Another source of errors in the experiment is the noise picked up 
by the receivers. Figure 4.9 shows the noise temperature, a measure of 
the amount of noise present, as a function of frequency. As can be seen, 
the frequencies used by the experiment, 2.225 KHz, 3.385 KHz, and 
5040 KHz, are in a portion of the spectrum with highly variable noise 
characteristics. The following is a calculation of the signal-to-noise 
ratio at the rocket for the 2.225 KHz system. 

The loss of signal due to the length of the propagating path for an 


isotropic radiator is 
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TABLE 4.4 


Receiver Output for 

Frequency Component 
dc 

2Kv 

r 

AKv 

r 

2nKv 

v 


- -1 and R = 0 


Altitude 


4 Q 

IT 


4£ 

3tt 


4g 

15tt 

4Q 

(4n 2 -l)ir 



AMPLITUDE (arb. units) 



FREQUENCY (Hz) 


Figure 4.8 Frequency spectrum of receiver output signal 
for Nike Apache 14.270, prepared by K. L. 
Miller [Edwards, 1973]. 







(4.5) 


L - 10 log 


4 tt xh 

\ 


under the assumption of no absorption. If the rocket is at 100 km, the 
loss is 50 dB at 2.225 KHz. The transmitter power (Py) is originally 
about 5 watts for the extraordinary signal which is 7 dBw. The combined 
gains of the antennas (G) are about 2 dB. The resulting signal power at 
the rocket is then 


P R - P T + G - L (4.6) 

which is -41 dB in this case. 

The total noise power received is 

N = K' + T + B (4.7) 

in which K' is 10 log k (Boltzmann's constant), equal to -229 dBw/°K Hz, 

T is 10 log T (noise temperature) and B is 10 log BW (bandwidth) . For 
a good payload receiver, BW is 1500 Hz. The temperatures chosen from 
Figure 4.10 are about lOO^K at minimum and 10^°K at maximum. These 
result in a maximum noise power of -100 dB and a minimum of -180 dB. 

Comparing the signal power to the noise power, the signal-to-noise 
ratios of 139 dB and 59 dB result. Clearly random atmospheric noise is 
not an important factor over the 30 dB range of absorption the experi- 
ment is designed to operate over. 

In addition to random atmospheric noise the signal may include 
unwanted radiofrequency transmissions. The receiver's bandpass 
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characteristic is poor in this respect. Examination of bandpass curves 
of receivers used in the payload shows a slow falloff and extraneous 
peaks. This can become important when the frequencies of the experi- 
ment are not far removed from those used by others. 2225 KHz is 
located in a frequency range used for radiotelegraph and ship-to-shore 
communications. Both of the upper frequencies used are within the 
tropical broadcast bands, used by countries near the equator because of 
poor propagation at normal broadcast band frequencies, as listed in 
Tables 4.5 and 4.6. Part of the hazard presented by these other users 
can be reduced by monitoring the receiver output prior to launch with 
the transmitters momentarily shut off. If there are any strong signals 
present, the rocket should not be launched unless poor performance of 
the differential absorption experiment is risked. 

Another area in which other users of the airways may present a 
problem occurs because the receivers, when placed in a payload, receive 
signals at frequencies other than those for which they are designed. 

Table 4.7 shows the location of peaks outside the passband and the 
amount of attenuation at these points. It must be pointed out that 
this is not a particularly bad receiver and that all units tested 
showed this type of response. The frequency locations of these peaks 
do vary between receivers. Several of these fall in the broadcast band 
and if there is a nearby station it could affect data. This problem 
will show itself, if the above test for other stations near the operation 
frequency is done. 
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Frequency 

3380 


3385 


3390 


TABLE 4.5 

Stations Near 3.385 MHz 
Power 

in Kilowatts Country 


18 

100 

1 

.25 

1 

.075 

.3 

10 

.6 

4 
1 
1 
1 

10 

120 

.2 

.25 

5 


Mali 

Malawi 

Guatamala 

Bolivia 

Peru 

Indonesia 
Indonesia 
Sri Lanka 
New Guinea 

Dominican Republic 

Brazil 

Peru 

Venezuela 

Afghanistan 

China 

Panama 

Bolivia 

Ecuador 


3 


Peru 
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TABLE 4.6 

Stations Near 5.040 MHz 


Frequency 

Power 

Country 

5037 

1 

Angola 


50 

USSR 


30 

Central African Rep. 


1 

Brazil 


1 

Columbia 


5 

Ecuador 


V .5, 1 

Peru 

5038 

30 

Central African Rep. 

5039 

20 

Sudan 

5040 

50 

USSR 


50 

Burma 


none given 

China 


.05 

Indonesia 


.15 

Haiti 


1 

Ecuador 


10, .5 

Peru 


1 

Venezuela 

5041 

10 

Guinea 


1 

Phillipines 


1 

Cook Islands 


5 

Bolivia 


1 

Brazil 


1 


Peru 
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TABLE 4.7 

Extra Signals Detected by 
Receiver #106 in 5040 Payload 
Frequency Attenuation 


556 

-26 dB 

585 

-52 dB 

601 

-52 dB 

635 

-52 dB 

674 

-52 dB 

695 

-52 dB 

741 

-30 dB 

856 

-46 dB 

927 

-47 dB 

1112 

-23 dB 

1236 

-47 dB 

1309 

-43 dB 

1483 

-34 dB 

1853 

-42 dB 

2780 

-38 dB 

3177 

-28 dB 

3705 

-54 dB 


4449 


-20 dB 
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5. RANGE AND ACCURACY OF DIFFERENTIAL ABSORPTION DATA 

In view of what has been discussed in the previous chapters about 
errors, it is appropriate to define a basis upon which the useful range 
of data may be determined. Chapter 4, Section 2, established an upper 
limit on the experiment, as a result of the polarization error. Here the 
lower limit of useful data will be fixed and some discussion of the 
accuracy of the measurements will be undertaken. For this purpose. Tables 

5.1 to 5.3 present a summary of the results of the rate recovery program 
for three flights. 

5.1 Lower Limit of Differential Absorption Data 

The problem is to estimate the lowest differential absorption rate 
which can be considered to be a significant measure of real absorption. 

To get a measure of the random errors, several seconds of data are 
processed before absorption starts . This is then used to compute the 
standard error of measurement. Table 5.4 shows the errors for the 
flights under discussion. The criteria for determining the starting 
point for good data is to select the first data point greater than the 
standard error after which there are two or more data points also 
greater than the standard error. This point up to where polarization 
error takes over is the range of useful data. 

Figures 5.1 to 5.3 show the probe current and only the useful dif- 
ferential data set by the procedures above and in Chapter 4. The 
polarization was not checked for the higher frequency in each experiment 
so the cutoff is taken as the first data point of value less than the 
standard error. 
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TABLE 5.1 

14.511 - 2.225 and 3.385 KHz 
Differential Absorption Rates 

Time 2.225 rate 3.385 rate 

after launch (dB/sec) (dB/sec) 

49 .07 .04 

50 .05 -.19 

51 .13 .27 

52 .25 .27 

53 .27 -.26 

54 .27 .36 

55 .33 .45 

56 .60 .04 

57 .96 -.35 

58 1.79 .89 

59 3.06 .92 

60 8.49 2.07 

61 2.39 2.75 

62 .22 2.25 

63 2.78 

64 2.56 

65 .73 

66 .91 

67 1.26 

68 1.12 

69 .62 

70 1.42 

71 .02 

72 1.31 

73 1.73 

74 2.10 

75 -1.16 

76 1.34 

77 .30 

78 -.46 

79 .62 

80 .11 


C 
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TABLE 5.2 

14.513 - 2.225 and 5.04 MHz 
Differential Absorption Rates 


Time 

2.225 MHz 

5.04 MH 

after launch 

(dB/sec) 

(dB/sec 

40 

.12 


41 

.10 


42 

.04 


43 

-.08 


44 

.02 


45 

.10 


.46 

-.04 


47 

-.04 


48 

.20 


49 

.02 


50 

-.12 

.18 

51 

.36 

.18 

52 

.44 

-.71 

53 

-.12 

.37 

54 

.58 

-.16 

55 

.74 

-.05 

56 

1.04 

.27 

57 

1.56 

.16 

58 

2.96 

.38 

59 

2.84 

.07 

60 

3.70 

-.05 

61 

3.21 

.49 

62 

2.53 

.02 

63 

2.10 

.18 

64 

1.73 

.22 

65 

.79 

-.02 

66 

.42 

-.05 

67 

.26 

.75 

68 

.07 

-.51 

69 

.07 

.42 

70 

.31 

.33 

71 

.52 

-.18 

72 

.57 

.58 

73 

.83 

.47 

74 

.87 

.37 

75 


.47 

76 


1.30 

77 


.84 

78 


.82 

79 


1.57 

80 


13.75 

81 


-.68 

82 


1.33 
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TABLE 5.2 

14.513 - 2.225 and 5.04 MHz (cont.) 

Time 2.225 MHz 5.04 MHz' 

after launch (dB/sec) (dB/sec) 

83 .12 

84 -.20 

85 -.27 

86 - . 88 

'87 4.56 

88 -1.25 

89 1.75 

90 -2.58 
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TABLE 5.3 


Time 


42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 


14.514 - 2.225 and 5.04 MHz 
Differential Absorption Rates 

2.225 MHz 5,04 MHz 

(dB/sec) (dB/sec) 


-.11 

.07 

.05 

.02 

-.16 

.10 

.13 

-.03 

.11 

0 

.03 

0 

.57 
. 16 
,70 
1.16 
1.25 
2.33 
3.09 
2.88 
2.68 
2.53 
2.22 
1.36 
1.41 
1.65 
1.44 
1.12 
-.03 
- .06 
.29 
-.09 
.92 
1.01 
.71 
0 
0 
0 


-.05 

.02 

.47 

-.•46 

-.20 

-.07 

.38 

.24 

.02 

-.15 

.07 

.46 

.11 

.36 

-.02 

-.11 

-.09 

.26 

.15 

0 

.22 

1.11 
-.31 
.58 
.26 
.51 
.64 
- .46 
.46 
25.58 
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TABLE 5.4 

Lower Limit of Data 




Standard 

First Significant Data 

Flight 

Frequency 

Error 

(dB/sec) 

Time from 
Launch (sec) 

Altitude 

(km) 

14.511 

2.225 

.07 

51 

61 


3.385 

.29 

58 

74 

14.513 

2.225 

.16 

54 

65 


5.04 

.33 

72 

91 

14.514 

2.225 

.09 

73 

65 


5.04 

.26 

73 

92 


v 
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Figure 5.2 Profiles of probe current (upper scale) and electron density 
(lower scale). 
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Figure 5.3 Profiles of probe current (upper scale) and electron density 
(lower scale) . 
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5.2 Error in Calculated Electron Densities 

The shaded area in Figure 5.4 shows the range over which the electron 
densities for 14.511 can fall if the errors of up to the standard error in 
Table 5.1 in the 2225 KHz system are allowed. As would be expected, as 
the rocket goes further and the rate of absorption increases, the random 
error represents a smaller portion of the total data. It will be shown 
below that standard errors previously calculated are not valid above the 
first few data points, but are a reasonable estimate. Figure 5.5 shows 
the standard error of each second of data as calculated from the rates of 
absorption taken one tenth of a second at a time. The magnitude of the 
error is not in agreement with that calculated in the last section. The 
reason for this is that the errors are not independent because of the 
feedback in the system. If the servo moves too fast, resulting in a 
high absorption rate being recorded, the feedback will result in some 
data where the rate is slower as the servo returns to its proper 
setting. The result is that the time over which data are averaged 
affects the statistics in that the longer the time the more closely the 
errors approximate independent random values. For this reason, the 
error values in Figure 5.5 are of larger magnitude than when whole 
seconds of data were taken at a time. To limit the effect of real ab- 
sorption, the errors are calculated from a straight line fit of the 
data generated for that second. Clearly the result shows that as the 
rocket penetrates the ionosphere the errors increase. The error is 
increasing at a rate near a factor of two for every ten kilometers of 
ionosphere traveled through. If this factor is used starting from the 
first valid point, the error range shown by the dotted line in Figure 
5.4 results. 



NIKE APACHE 14.511 
C = 3.88 X I0 5 p 



J l I 1 I I I J l 1 L 

I0 2 

ELECTRON DENSITY (cm' 2 ) 


Figure 5.4 Estimated error in electron density. 
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5.3 Collision Frequency 

Mechtly [1974] discussed the problem of analyzing differential ab- 
sorption data in the absence of Faraday rotation. It will be recalled 
that under these conditions a collision frequency model of the form 
V m = C x 10^p is used. Figure 5.6 shows the results of using C - 3.88 
and C = 7.3 on the data from 14.511. Looking at Figure 2.1 one notices 
that using the Sen-Wyller formulation, the absorption is nearly inde- 
pendent of collision frequency at 72 km. This is confirmed by the 
crossing of the two curves in Figure 5.6 near this altitude. 

At points where the data starts, near 62 km, the error between the 
curves resulting from the choice of C, is greater than that due to 
inaccuracy in the measurement of the absorption rates.. One point that 
should be noted is that at the 72 km point both errors in measurement 
and of collision frequency selection are only a small fraction of the 
total electron concentration. This means that if the probe experiments 
electron concentration to current ratio is constant in the lower iono- 
sphere the 72 km point becomes an excellent calibration point. 



NIKE APACHE 14.51 
0 C* 3.88X10* p 

© C = 7.3 X I0 5 p 


ELECTRON DENSITY (cnrf 3 ) 


Figure 5.6 Electron densities using different collision frequency models 
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6. CONCLUSIONS AND RECOMMENDATIONS 

6.1 Areas of Experiments Usefulness 

As a result of careful examination of experimental equipment and pro- 
cedures, it is felt that the objective of improving the experiment to a 

-3 

point where electron concentrations on the order of 10 cm are measured, 
is not obtainable, without major system changes. The level to which the 
experiment can be used is a function of the absorption the wave goes 
through. 14.511 which was launched on a day of high absorption allowed 

_3 

the experiment to be used down to about 100 cm . The other day shots 

examined, 14.513 and 14.514, for which the absorption was not as high, 

-3 

were useful to levels of 200 cm . For nighttime shots with little 
absorption, the experiment is generally rendered useless by reflections. 

As a result of the analysis of where the data can be used, the 
2.225 MHz system is found to be useful in the day generally only between 
65 and 75 kilometers. The use of the 5.04 MHz system is useful only over 
an area of about five kilometers. The 3.385 MHz, not generally used in 
the day, is found to work reasonably well giving much more information 
than available from the 5.04 MHz system. 

6.2 System Change Recommendations 

The following are suggestions for changes in the system which will 
lead to improved operation of the differential absorption experiment. 

None of these changes will degrade the operation of the Faraday rotation 
experiment. 

6.2,1 Polarisation. To stop the servo system from following the 26 
signal generated by the error signals, filters can be employed. Two 
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possibilities exist.- A notch filter set to eliminate any signal at 26 can 
be used, or a bandpass filter centered at 26 + 2u>^ with enough attenuation 
at 26 to effectively make it 0. For either system, the filter would have 
to be bypassed for ground tests and during early stages of flight before 
the spin rate becomes appreciable. The bandpass filter has a disadvantage 
in that different characteristics would be needed for rockets using dif- 
ferent spin rates. The filter must be of very narrow characteristics in a 
slow spin rocket since it would be desirable to have 40 dB of attenuation 
in, the 26 signal only 2 u radians/sec away from the center of the pass- 
band. Some of the bandpass filter’s faults could be overcome by the use 
of a tracking filter which will follow the signal as the spin rate 
changes. These could be set at 26 at launch and would then follow the 
signal as the spin rate builds up, providing constant system lock. 

For the notch filter only one need be designed centered at 26. It 
should provide 20 dB of attenuation at 26 and essentially none at 26 + 2 
when the smallest spin rate is present in the area of data considered. 

In addition to filtering to separate the 26 and 26 + 2^ signals, an 
improved filter to separate the dc component is needed. The servo has 
been noted to follow the 2w^ modulation present. The present filter is a 
first order filter at 1 Hz. It would be desirable to replace this by one 
with an upper 3 dB point at 2 Hz and with attenuation of. 40 dB at about 
8 Hz. This dc filter would most likely improve the standard error of 
the measurements. 

It may be advisable to change the value of 6, which is now 250 Hz, 
if a filter is used. By changing 6 down to 50 the $'s of the filter 
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sections can be reduced by a factor of 5 making stability of the filter 
less of a problem. This would require, however, redesign of the modula- 
tion detector. ■ 

Another approach to the problem is to improve the actual polarization 
characteristics. This is not deemed practical, as an improvement of a 
factor of about 50 would be needed. This type of improvement could not 
be obtained with a simple antenna system required for field operation. 

6.2.2 Receivers . -■ As noted in Chapter 4 the receivers have many 
undesirable characteristics. These characteristics have been worse in 
the newer receivers. Since the receivers were originally designed in 
the early 60’ s many improvements in techniques have been made. It is 
suggested that new receivers be designed. These new receivers need not 
be as sensitive as those presently used. Any sensitivity beyond that 
necessary to detect the propagating wave, after going through 30 dB of 
attenuation, will only result in increased noise pickup. 

6.3 Data Reduction 

Most of the improvements in this area are under way. The data now 
comes digitized on tapes, eliminating that step in the processing. What 
needs to be done now is to streamline the processing by reducing the 
amount of human data handling. A start in this area is being made by 
having each program punch its results onto cards with the appropriate 
data for the next step. The goal should be to design the processing so 
that two sets of cards are generated, one with differential absorption 
data and trajectory, and the other with Faraday rotation data and tra- 
jectory. These could then be used in a single program to generate elec- 
tron densities. 
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6.4 Further Test Suggestions 

The effect the telemetry system has on the accuracy of the experiment 
was not studied. It would be desirable to make a few tests while the 
rocket is on the ground to measure the accuracy of the telemetry system. 
One test would be to operate all systems as if under flight conditions 
and record the servo systems movements. This will test the noise gener- 
ated by the transmitter to rocket to telemetry station to van loop. 

Next simulate the expected receiver output at the VCO’s in the payload 
to eliminate the transmitter to rocket section of the loop. Finally, 
simulate the expected inputs to the data lines at the telemetry station 
to test the line and the modulation measuring equipment in the van. 

If the test of the loop shows errors in recording on the order of 
.07 dB/sec when analyzed as real data, the loop may be the source of 
most of the experiment’s errors. The further tests will help evaluate 
where the most improvement can be made. 
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